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Abstract 

For over 25 years it is known that the roll structure of electroconvection (EC) in 
the dielectric regime in planarly aligned nematic liquid crystals has, after a transition 
to defect chaos, the tendency to form chevron structures. We show, with the help 
of a coarse-grained model, that this effect can generally be expected for systems 
with spontaneously broken isotropy, that is lifted by a small external perturbation. 
The linearized model as well as a nonlinear extension are compared to simulations 
of a system of coupled amplitude equations which generate chevrons out of defect 
chaos. The mechanism of chevron formation is similar to the development of Turing 
patterns in reaction diffusion systems. 


1 Introduction 


Understanding defect chaotic states of pattern forming systems is presently 
one of the important goals of research in pattern formation |^. Important 
questions relate to the physical quantities characterizing properly such a state 
[^|-^ and to the possibility to observe and understand transitions between 
different types of defect chaos @,0]. 

An interesting candidate for such a transition is the formation of chevron 
patterns occurring in particular in the dielectric regime of electroconvection 
(EC) in thin layers of planarly aligned nematic liquid crystals In¬ 

creasing the ac voltage applied across the layer (above the ’’cutoff frequency”) 
one hrst observes a (nearly) periodic pattern of (narrow) convection rolls with 
wavevector parallel to the director orientation imposed by surface treatment 
of the plates forming the upper and lower boundaries of the layer. Then de¬ 
fects (dislocations) in the roll pattern start to appear, their density increasing 
as the voltage rises. Initially they move irregularly and their distribution is 
homogeneous. Defects carry a topological charge which is -|-1 or — 1 depending 
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on whether a roll “ends” or “begins” at that point. After a further increase 
of the voltage one observes the transition to chevron patterns, with defects of 
equal charge ordering along periodically arranged domains oriented parallel 
to the original roll direction. Simultaneously with the onset of the ordering 
process the increase of the number of defects with voltage becomes noticeably 


steeper 


Between the chains the convection rolls are rotated away from 


their original direction alternatingly to the left and to the right, which gives 
the whole structure a “herringbone” like appearance. 

In this work it will be argued, that the tendency of a 2-d anisotropic pattern¬ 
forming system to form chevron patterns is a rather general feature depending 
in particular on symmetry properties. As a matter of fact chevron patterns as 
shown in Fig. (for details see Sec. can easily be observed in simulations of 
the recently proposed generic amplitude equations (|^,^) |^. This description 
was inspired by EC in the low-frequency conductive regime in homeotropically 
aligned nematics with negative dielectric anisotropy where one starts from a 
situation without external anisotropy (director in the z direction perpendicular 
to the slab). Increasing the voltage one hrst has a spatially homogeneous 
Freedericksz transition where the director bends away from the z direction 
thereby singling out spontaneously a direction c in the x-y plane. The isotropy 
may also be broken externally by applying a weak magnetic held parallel to 
the plain, c may then be expressed by the angle (p that is enclosed by c and 
the held. At higher voltages there is an instability to EC with a critical mode 
corresponding to rolls with wavevector parallel to the in-plane director c (in 
the simplest case). A weakly nonlinear description then has to incorporate 
the roll mode as well as the undamped mode corresponding to (inhnitesimal) 
rotations of c (the Goldstone mode). Equations (|l|,|D below describe such a 
situation. In fact chevrons have recently been observed in homeotropic EC 
|T^ . Other systems where chevrons were observed are dielectric EC under 


oblique boundary conditions H. EC in nematic polymers and sometimes 


also in the conductive regime of EC in planarly aligned ordinary nematics |19 


After introducing our notation and dehnitions in Sec. we develop a linear and 
(weakly) nonlinear model of chevron structures in Secs. § and In Sec. ^ some 
predictions of the model are tested quantitatively against numerical simula¬ 
tions of Eqs.(mD. Section ^ discusses experimental and theoretical questions 
in a more general framework. 


2 Formal Setting 


The derivation from symmetry arguments of the coupled amplitude equations 
for systems with spontaneously (and almost spontaneously) broken isotropy is 
presented in . (There the coefficients are also calculated from hydrodynam- 
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ics for EC in homeotropically oriented nematics.) The system should bifurcate 
supercritically to a roll pattern with wavevector parallel to c = (cos (p, sin (p) 
(normal rolls) when some external control parameter e changes sign from neg¬ 
ative to positive. For small e and (p and after some rescaling the equations 
take the form 


rdtA ={l + dl + {dy - iipf + ipyip^y - \ A\^)A, (1) 

d,^={K,dl + dl-H^/e)^+ ( 2 ) 

T{—iA*{dy — iip)A + C.C.). 

In Eq.®) the derivative operator dy operates only on A and ip^y := dip/dy. 
All coefficients are real, r and are positive. A is the complex amplitude 
of the patterning mode. A = const, corresponds to the most unstable linear 
mode at = 0. Similar to the usual real Ginzburg-Landau equation a factor 
~ has been taken out of the physical length scales, a factor ~ out 
of time, and a factor ~ out of A and (p. The spatial extensions of the 
plane L^, Ly should be large {Lx,Ly 3> 1). Whenever needed we will impose 
the convenient periodic boundary conditions. The coefficient containing 
describes the (small) external perturbation of isotropy. (We use this notation 
since H is typically a (scaled) magnetic field.) For sufficiently large (or, 
equivalently, small £ while H ^ 0) one has a stable band of stationary, spatially 
periodic solutions of (|^,^. All such solutions are unstable for = 0 if T is 
negative, which appears to be typical for nematics. In simulations one typically 
observes a transition to a defect chaotic state when h? := H‘^/{—2Te) drops 
below a critical value of 0(1) |1l5| . 

Writing the complex amplitude A as A = |A| e*® defines the phase modulations 
6 of the underlying stripe pattern up to multiples of 27r. This degeneracy plays 
no roll as long as only 

ve = lm|^| (3) 

is considered (and A 7 ^ 0). V 6 * is the deviation of the local wavevector from the 
most unstable wavevector at (p = 0 (up to the rescaling done in Eq. (|^)). The 
wavevector of the roll pattern and the director c are parallel if P := dyO = ip. 

A topological defect in the stripe pattern corresponds to a simple zero of the 
complex amplitude A. Its topological charge is / dr*- V6'/27r, where the path 
of integration is a small loop encircling the defect in the positive sense. 

We define defect densities 


n±{r) = ^5(f-f±j), 
j 


( 4 ) 
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where the sum is over all positively or negatively charged defects at the po¬ 
sitions f±j, respectively, and S{-) is the Dirac 6 function. The total defect 
density is n := n+ -|- n_ and the topological charge density is p := — n_. 

One has 

f V6{f) df = 271 ( p(r) dO (5) 

an n 


for any area D. 

From an abstract point of view one can define chevrons as a periodic modu¬ 


lation of p, with a wavevector parallel to the x axis [|r^. Hence we will only 


look at modulations of the defect chaotic state in the x direction and average 
all equations along ?/, which will be indicated by an overbar “ := L~^ f ■ dy. In 
particular one hnds a topological condition 


d^P = 27rp, 


( 6 ) 


by differentiating both sides of Eq. (i for a rectangular area D = (0...^) x 
(O..Ly) with respect to Obviously P can only change through defect mo¬ 
tion. In fact, for fixed and (locally) constant p, where Eq. m reduces to the 
simple Ginzburg-Landau equation, defects always move such that the (local) 
wavenumber approaches the value of highest growth rate 
the condition P = p. 


211 , which includes 


3 Linear Model 


A simple model for the coarse-grained dynamics of the defect-chaotic state 
helps to understand chevron formation as a linear modulation instability of 
homogeneous defect chaos. Led by symmetry considerations |]^ we propose 
a model equation for the conservation of topological charge 


dtp + -DdxP + a{P -p) 


= 0 . 


( 7 ) 


The diffusion coefficient D and the “conductivity” a are phenomenological 
constants. We do not calculate them here and assume both to be positive. This 
means that diffusion of defects acts to minimize topological charge imbalance 
and systematic drift acts to reduce \P — p\. To get a corresponding equation 
for p we average Eq. (|^) over y. 

dtip = (K,al + 2rh^)<p + 2T\A\lg(P - Tp). (8) 
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The form of the last term follows from isotropy for small P — Tp^ is a 

suitable constant. Equations (§|,0)0) describe dynamics for small deviations 
from p = 0 on scales much larger than 

Using this model one can calculate the stability of homogeneous defect-chaotic 
solutions against periodic modulations of P and ^ in a: direction. After apply¬ 
ing dx on Eq. (|[) and eliminating P through (||), Eqs. take the form of 
a linearized reaction diffusion system for p and dx^- Accordingly one gets a 
homogeneous Hopf bifurcation or a steady state, spatially periodic (diffusive) 
Turing instability or a simple non oscillatory, homogeneous bifurcation as the 
Erst instability (see e.g. ). 

For ATs/D > 1 and 1 < -T\A\ls/7ra, or K^^/D < 1 and {D/KsY/^+iKs/Dy/^ < 

2(—riAIgg/vrcr)^/^ one has the Hopf bifurcation at -|- ira/T. The 

Hopf frequency is = 47r(T(—r|A|gg — vra). 

If {D/KY^!^ + > 2(-r|A|2g/7ra)V2 and K^/D < -^A^^/^Ka the 

chevron (“Turing”) pattern with critical wavenumber k1 = {—A'KT\A\‘lQa/DKY)^/‘^ — 
2'Ka/D appears first at = (|A|efr — {7rKsay^‘^{—DT) ^/^)^ (see also Fig. |I]). 

None of the two instabilities occurs in the remaining case that K^/D > 
-mi^/ TTCT and 1 > —r|A|gg/7rcT. At = 0 the mode of homogeneous 
rotation (Goldstone mode) invokes instability. 

According to this model a crucial effect is the diffusive contribution in Eq. (^. 

It invokes a wavevector mismatch [P — ip) for modulated P and p which would 
otherwise be leveled out by defect motion. Due to this mismatch the repul¬ 
sive “torque” on p by the roll pattern (described by the term containing T in 
Eq. (^) can then drive chevron formation. It should be noted that no particu¬ 
lar assumption about the interaction of single defects was made. By inspection 
one finds that the model is robust against additional terms in Eq. (|^, as long 
as modulations of P are weakened compared to modulations of ip. For ex¬ 
ample, a consistent truncation in powers of k should include a fourth-order 
diffusive term in Eq. (|^) that has been avoided for the sake of simplicity. 


4 Nonlinear Model 


In the final, saturated chevron pattern the defect density n becomes strongly 
modulated with half the period of the chevron pattern. The defects accumulate 
along “chains”. This can be understood with ideas similar to those describing 
a p-n junction in a semiconductor. In the regions along the chains a high 
density of positively or negatively charged defects is enforced by the strong 
bent {dxp) of the director. This corresponds to doping the p and n regions of 
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the semiconductor with charge carriers. The density of the oppositely charged 
species is strongly suppressed due to high recombination probability. As in the 
depletion layer of the p-n junction, in the region between the chains the total 
density of defects n is reduced compared to the “chain” region. 

The coupling between p and n may also provide an important contribution 
to the nonlinear saturation of the chevron mode. Assume that defects are 
homogeneously created at a constant rate uo/r„, while they are annihilated at 
a rate proportional to n+n_, such that for the homogeneous defect chaos one 
has a time-averaged defect density (n)^ = 2no (In this section we only consider 
quantities averaged over y and drop the overbars.) Ignoring correlations in the 
nonlinear terms we obtain 


dx3± + dtn^ = - {nl - u+u_), 


TnUo 


(9) 


where j± are the defect currents. Since = Uq — n^/4 -|- p^/4 an 

excitation of the chevron mode increases the equilibrium number of defects. 
One should expect that with increased defect density n the conductivity a 
becomes higher (e.g. a = mn with some mobility constant m). This weakens 
the relative importance of the diffusive term in Eq. ( 0 ) so that the chevron 
mode saturates. In order to get a quantitative estimate of this effect we calcu¬ 
lated the weakly nonlinear m steady state solution of Eqs. ii) , assuming 
currents of the form 


j± = -Dd^n± ± m n±(P - (p), 


( 10 ) 


and allowing for a modulated value of |A|gg by substituting 


l^le\ - l^le\ (1 - on - - C,{P - 


( 11 ) 


in Eq. (||). 

For periodic solutions with wavevector k one obtains, up to a phaseshift. 



(12) 


(13) 

n-2nQ+ ^ ^ cos^fcx-F 0(h^ 

loTT^TT-ow 

(14) 


with 
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( 15 ) 


G = —Dk^lAllg {64:Ti‘^rn^nl{Dk‘^ + 4:TrmnoY)~^ ^ 

[27rm^no + 

ci{Dk^m^no + ATim^nl) + 

C2{4:Dk‘^m‘^nl + IGnm^nl) + 
c^iyinD’^k'^mriQ + ). 

Here = h‘^{k) is the value of where modulations with wavevector k 
become unstable. 


Note the following points: The result is independent of r„. However, this de¬ 
generacy can be removed, e.g. by introducing cross diffusion terms —D 2 dxn^ 
in Eq. (W). Another consequence of the particular form of Eqs. (^,]T) is the 
simple relation 


^2fc,0 


\pk,o? 

4no 


(16) 


between the Fourier modes of p and n [we define the Fourier transform / of 
a function / by f{x,y) = Y.k,i fk,i^^pi{kx + ly)]. To lowest order one has 
n — 2no ~ cos^ kx^ so that right between the chains there is no net effect of 
the pattern on the value of n. Both results should allow a simple comparison 
with experimental chevrons. 


5 Comparison with Simulations 


As pointed out in the introduction, chevron patterns are easily generated by 
simulating our “microscopic equations” (|l],|D. Figures and ^ show snapshots 
of the steady state for different sets of parameters. 

For the semi-quantitative test of our model we used a pseudo-spectral algo¬ 
rithm on a 128 x 128 grid with spatial extensions = 159.63, Ly = 60.125 
and periodic boundary conditions. 

The set of parameters in (|I],|[) was chosen in such a way that the chevron 
bifurcation sets in with a small wavevector. Moreover we made sure, that the 
oscillatory mode is absent. These conditions are satished by keeping 

T = 1.53558, I3y = 1.07013, 

^ (17) 

As = 0.180594, T = -0.0304092 


fixed, while taking as a control parameter. We remind that the chevron 
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amplitude increases by lowering h^. The condition of small implies also 
h‘^{kc) to be small and consequently the supercritical range is rather small, 
too. 

The most basic assumption of the model is that, to linear order, the dy¬ 
namics of P and ip may be isolated from the remaining degrees of free¬ 
dom and can be described by a set of linear PDEs, even though the full 
dynamics of A is strongly nonlinear. The coupling to the other degrees of 
freedom should be describable by adding noise. Close to the chevron insta¬ 
bility, where one of the branches of eigenvectors of these linear PDEs is 
strongly excited by the noise, there should be a strong correlation between 
the Fourier modes of P and Figure |] shows the correlation coefficient 
r{k) = Re{Pk,o(pl^o)t/i{\Pk,o\'^)t{\^k,o\'^)ty^^ for = 0.025. Clearly there is 
essentially perfect correlation for k < 0.24. The sharp drop of r{k) at /c ~ 0.3 
may be partly due to a change of the relative sign of Pk^o and pk,o in the linear 
mode. 

In order to test the usefulness of the particular form (|^||) of the PDEs we 
hrst developed methods to “measure” a and D in simulations. To hnd a, the 
Fourier mode po^o of p is hxed at some hnite value po, while all other modes 
evolve according to Eq. (||). This was implemented by resetting <^o,o to po after 
each time step. The simulation is run until a steady state is reached. Then 
Po is set equal to zero. The relaxation rate of the global average of P equals 
27rcr. Figure ^ shows the average value of P as a function of the time At after 
switching (po to 0. Within the accuracy obtained, the decay is exponential. 

Similarly, in order to measure D, a single pair of long wavelength Fourier 
modes pko,o, of p was held hxed at a hnite value po. D/a can be 

calculated from the steady state average value of Pfco,o by setting dtp = 0 in 
Eq. (|^. The values for a and D obtained for diherent values of po and ko 
were consistent. As a simple estimate for |A|gg and no we took the spatial and 
temporal average of |Ap and n/2. 

The numerically calculated coefficients at = 0.025 are 
|AL\ = 0.7996(4), no = 0.00509(10), 

\ -Lo J 

D/a = 38.0(1.3), a = 0.01330(73). 


Error estimates contain only stochastic contributions. From the linear model 
one then calculates the threshold h/. = 0.037(5) and the critical wavenumber 
kc = 0.214(5). 

Critical slowing down and strong noise make a precise determination of the 
threshold in simulations difficult. Since we expect to be close to zero, where 


the Goldstone mode becomes unstable, the straight forward method of extrap¬ 
olating the supercritical modulation amplitude to zero cannot be applied. In¬ 
stead we looked at the subcritical excitation of linear modes by noise. Without 
nonlinear interaction one expects a relation of the type (|pfc,oP)j ^ -Sk{h‘^), 

where Sk{h?) is the growth rate of the linear mode at wavevector k. Figure ^ 
shows some numerical results for (|pfc,oP)t at various values of k. The dot¬ 
ted line is a linear £t to fc = 0.19680, which implies a threshold value of 
^cnum ~ 0 . 01 . At however, the k = 0.19680 mode itself is already 

strongly suppressed by nonlinear interaction with other modes. We believe 
that the discrepancy between he and hc^num is to a large part an effect of the 
combination of nonlinearity and strong noise. For further discussion see next 
section. Figure shows the fully excited chevron pattern at = 0.005. 

We tested two aspects of the nonlinear model. Firstly, we hnd Eq. (pT]) to be 
well satished, as long as only one mode is active. In Fig. the distribution of 
pairs (|h 2 fc,o|) IPfc.oP) and the theoretical line for k = 0.19680 and h? = 0.025 
are shown. While \pkfi\^ is large the other modes are suppressed nonlinearly 
and the single mode approximation made in Eq. (PT| ) is legitimate whereas 
when |pfc,oP is small the influence of competing modes becomes noticeable. 


Secondly, we tested the prediction of Eq. (B) with (p!5[) and (|^) for the chevron 
amplitude. Due to the intricate situation at threshold our simple result cannot 
give more than an order of magnitude estimate. The solid line plotted in Fig. |^, 
shows the prediction of formula (^. It was calculated using coefficients as for 
the linear model with k = kc, m = a/2no = 1.31(8) and ci, C 2 , C 3 = 0 without 
any htting. It is easily seen that in our case the terms containing ci, C 2 , C 3 give 
only small corrections as long as Ci, C 2 , C 3 = 0 ( 1 ). 


The Hopf bifurcation predicted by our model can actually be found in simula¬ 
tions of Eqs. (1^,0. For example, at r = 2.19369, Py = 1.07013, = 0.338614, 

r = -0.0608184 and = 0.1 one gets \A\l^ = 0.7409(5), a = 0.01075(36) 
and D/a = 44.4(2.2). From this the Hopf threshold is found to be = 
0.186(19). The frequency of the fc = 0 eigenmode is o;^ = — 47 rF|A|gg(j — 
[r(|AL\ — h?) — -Kap = (0.0282(7))^ which is in good agreement with the 
maximum of the Fourier transform of the simulated time series of <^ 0,0 shown 
in Fig. 1^. 


6 Discussion 


As mentioned in the introduction the most striking case of experimental 
chevron formation occurs in the dielectric regime of electroconvection in pla- 
narly oriented nematics. At hrst glance one would not expect this system to 
fall into the class described by the theory, Eqs. (|^0), because here the horizon- 
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tal orientation of the nematic director n, which would represent the isotropy 
breaking degree of freedom, is hxed by the boundaries which have been pre¬ 
pared to align the director in one direction. However, in the dielectric regime 
the wavelength of the convection pattern is essentially independent of, and 
usually much smaller than the thickness of the layer d . Then the bound¬ 
aries can be reduced to a perturbative effect that may be described by the 
term H'^ip in Eq. (0). 


It is not easy to verify the validity of our model for planar dielectric EC. 
The most important prediction is the relation between the local wavevector 
and the director. A straight forward polarimetric measurement of the direc¬ 
tor making use of the birefringence of the material, however, is not possible, 
since the polarization axis of light passing the probe follows the director es¬ 
sentially adiabatically. The polarization axis is therefore determined by the 
boundaries. Some evidence for our model can be taken from the fact, that 
in the regime of strongly developed chevrons, where the local wavevector is 
turned by nearly 7r/2, new structures form inside the nematic, which super¬ 
impose with the chevron pattern. They allow an interpretation as disclination 
lines resulting from the fact that in the midplane the degenerated directions 
h and —n reconnect, which cannot occur at the boundaries El- 


Defects in experimental chevrons can align much better along chains than has 
been found in the simulations. This is presumably because higher-order terms 
and non-adiabatic effects that couple the defect cores to the underlying roll 
structure are not included in the amplitude equations. When experimental 
chevron patterns are strongly excited it is usually observed, that defects start 
to move along the chains, alternatingly ‘up’ and ‘down’, i.e. there is a corre¬ 
lation between dx^p and the y component of the velocity of a defect (ignoring 
its charge). This can not be observed in simulations of Eqs. (|l],^. In fact, it 
is forbidden by the accidental symmetry 93 — —(p, A ^ A* oi Eqs. (|^^. 
Inclusion of higher order terms would break this symmetry. In contrast the 
invariance under (p —>■ —(p, y —y, which follows from inversion symmetry 
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remains. 


We wish to emphasize that chevrons cannot appear directly from the homo¬ 
geneous (basic) state via a stationary supercritical bifurcation. Although this 
is consistent with most experiments there is some evidence that at sufficiently 
high frequencies chevrons can be observed directly at onset of the dielectric 
instability . This would mean that either the primary bifurcation becomes 
subcritical, which is not expected from theory, or there is a (“hidden”) bifurca¬ 
tion leading to an inhomogeneous state already below the dielectric instability. 
According to some older measurements the possibility of a hrst transition to 
a state with “wide domains” (with width ~ d) appears to exist El> 

but this phenomenon has not been cleared up. In any case this cannot be the 
general explanation. 
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It has been speculated that chevrons can be understood as a kind of inter¬ 
ference effect between the two most unstable dielectric modes, which differ in 
their z —> —z symmetry. Because of the character of the dielectric rolls (wave¬ 
length much smaller than d) the growth rates of these modes differ very little 
(as is the case for the higher z modes) [|^. The model has not been worked 
out, and can, in our view, not explain chevrons. 


It might be interesting to approach the chevron transition from the point of 
view of phase transition theory. We are dealing with a breaking of a continuous 
symmetry in two dimensions of the same symmetry class as the x-y model. 
Therefore the transition is, at least when nonadiabatic effects are neglected, 
possibly of the Kosterlitz-Thouless type . The model given here represents 
a simplified, mean-field type description. Consequently it is not surprising that 
the transition is delayed and difficult to pin down in the full “microscopic” 
theory described by Eqs. 


Finally we wish to point out that chevron-like structures have also been found 
in simulations of the anisotropic version of the well-known complex Ginzburg- 


Landau equation . The mechanism operative here is presently under inves¬ 
tigation. 
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Fig. 1. Regions in parameter space where chevron instability and Hopf bifurcation 
respectively are the first to occur when h? is decreased. The circle and the square 
correspond to the two parameter sets analyzed in Section 
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Fig. 2. Snapshot of a simulation of Eqs. (||,§) with r = 4.38738, (5y = 1.07013, 
K^, = 0.0564357, F = —0.0304092 in a rectangle of size 159.631 x 120.25. The upper 
image shows 1^41 encoded in a grey scale, the lower one shows the positions and 
polarity of defects as detected by the algorithm used throughout this work. 
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Fig. 3. Correlation between Fourier modes (see text). 
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Fig. 4. Relaxation of Po,o after a jump in (^o,o from (/?o = 0.1 to 0, averaged over 59 
runs. The dashed line is a least square fit. 



Fig. 5. Test of the linear and nonlinear model, see text. 
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Fig. 6. Samples of pairs (|n 2 fc,ol) l/Ofc,oP) at A; = 0.19680 and h? = 0.025 from numer¬ 
ical simulations and the relation given by (|l^) (solid line). 



Fig. 7. Snapshot of the steady state corresponding to h? = 0.005 in Fig. p. Top: 
values of (p (—2.33... 2.30) encoded on a gray scale; center: |A|; bottom: defect 
positions and polarities. 
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Fig. 8. Fourier transform of a simulation of Eqs. (|T|,§) in the oscillatory regime and 
the frequency of the most unstable mode of Eqs. (^-|^ (dashed line). The parameters 
are given in the text. 
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